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Abstract 

This paper gives an analysis of the periodic solutions of a ring of n os- 
cillators coupled to their neighbors. We prove the bifurcation of branches 
of such solutions from a relative equilibrium, and we study their symme- 
tries. We give complete results for a cubic Schrodinger potential and for a 
saturable potential and for intervals of the amplitude of the equilibrium. 
The tools for the analysis are the orthogonal degree and representation 
of groups. The bifurcation of relative equilibria was given in a previous 
paper. 

1 Introduction 

Consider a lattice of n nonlinear oscillators, with a periodicity condition or a 
ring of n oscillators. The lattice may come from the discretization of a non- 
linear PDE, such as a nonlinear Schrodinger equation, with different types of 
potentials. These equations have been used very often as models of several 
phenomena of physics of waves, in particular in nonlinear optics, see [2]. 

In previous papers, we have studied the problems of bifurcation of relative 
equilibria (see [I]) and bifurcation of periodic solutions for rings of masses ([5]), 
vortices and filaments ([5]). 

Although these problems are different as models of physical phenomena, they 
have in common several mathematical aspects, such as a similar linearization 
and the same type of symmetries. 

A relative equilibrium is a stationary solution of the system in a constant 
rotating frame. In particular, we consider one where all oscillators have the 
same amplitude but different phases. When these equilibria are localized they 
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are called breathers. In case of a periodic localized solution, one has a quasi- 
periodic breather, see [2] or [9]. 

Numerical studies of the bifurcation of relative equilibria for a circular lattice 
and a defect in one of the oscillators, for n = 7, is given in [12]. For the case of 
a Dirichlet condition, instead of the periodic problem, one has a very complete 
numerical analysis of bifurcation of breathers in . 

The linearization of the system at a critical point is a 2n x 2n matrix, which 
is non invertible, due to the rotational symmetry. These facts imply that the 
study of the spectrum of the linearization is not an easy task and that the 
classical bifurcation results for periodic solutions may not be applied directly. 
However, we shall use the change of variables proved in our previous paper, 
@], in order to give not only this spectrum but also the consequences for the 
symmetries of the solutions. 

The present paper is a continuation of [1]. Thus, we shall use the results in 
that paper, but we shall recall all the important notions. In our parallel papers, 
[5] and [BJ, we study similar problems for point masses and vortices. Although 
there are many similarities, in particular in the change of variables, the results 
are of a quite different nature. 

The next section is devoted to the mathematical setting of the problem, 
with the symmetries involved. Then, we give, in the following two sections, 
the preliminary results needed in order to apply the orthogonal degree theory 
developed in [5], that is the global Liapunov-Schmidt reduction, the study of 
the irreducible representations, with the change of variables of [4], and the 
symmetries associated to these representations. In the next section, we prove 
our bifurcation results and, in the following section, we give the analysis of the 
spectrum, with the complete results on the type of solutions which bifurcate 
from the relative equilibrium. 

2 Setting the problem 

Let us denote by qj(t) <E C the j'th oscillator, for j E {l,...,n}. The dNLS 
equations are 

iq 3 = h(\\qj\\ 2 )qj + (q j+ i - 2 qj + q^t), 

where h is the nonlinear potential. We wish to study a finite circular lattice, that 
is, a lattice of oscillators for j e {1, ...,n}, with periodic conditions qj = qj+ n . 

The solutions of the form qj = e uti Uj, with Uj constant, are called relative 
equilibria. In order to obtain the amplitude as a parameter, we need to change 
coordinates, with qj — fie Uj. In this manner, we have that the values Uj form 
a relative equilibrium when 

— Luuj — h(\/j,Uj\ 2 )uj + (uj+i — 2uj + Uj-i). 

Remark 1 Given that the lattice is integrable for n = 1 and n = 2, we shall 
look for bifurcation for n > 3. 
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The starting point is a relative equilibrium which looks like a rotating wave. 
We give a condition which needs to be satisfied by the potential for the existence 
of this rotating wave. 

Proposition 2 Define a,j = e 1 ^ , with £ = 2ir/n, then a — (a%, ...,a n ) is a 
relative equilibrium if 

LU = 4:Sm 2 (C/2)-h(p 2 ). 
The proof is simple and given in [4] . 

For the study of periodic solutions, we need to change to rotating coordi- 
nates: <7j = fie utl Uj, and the equation becomes 

iuj — ujUj = h(\/j,Uj\ 2 )uj + (uj+i — 2uj + 

Changing to real coordinates, one has, for uj £ M 2 

Jiij = LOUj + h(\fiUj\ 2 )uj + (uj + i — 2uj + Uj-i), 

where J is the standard symplectic matrix. 

Set u = (ui, •••) u n ) T be the vector of positions and let J — diag(J, J), 
then the vectorial form of the equation is 

Jii = W(it) with 
1 n 

i=i 

where H(u) is the function which satisfies V-ff (it) = uju + h(\fj,u\ 2 )u. 

Since we are looking for periodic solutions, let us define x(t) = u{t/v). Then, 
2tt/v solutions for it are 2tt solutions for x of the equation 

f(x) = -vjx + VV(x) = 0. 

The operator / is defined from H^(R 2n ) into ^(R 2 "). 

Definition 3 Let S n be the group of permutations of {1, n}. One defines the 
action of S n in M. 2n as 

P(l)(xi,--,x n ) = (x l{1) ,...,x l{n) ), 

and the action of 8 G 50(2) as 

p(9)x = e" J6 x. 

Let Z„ be the subgroup of permutations generated by C,(j) = j + 1 modulus 
n. The gradient \7V is Z„ x S'0(2)-equivariant, that is it commutes with the 
action of the group, and the map / is T x S ^equi variant with the abelian group 

r = Z„ x SO(2), 
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where the action of S 1 is by time translation. 

Now, the infinitesimal generators of S 1 and T are 

d d 
Ax = — \u, = Qx(t + ip) = x and A : x = -7rr\e^~ J6 x = -Jx. 
dip aO 

Since V is T-invariant, then the gradient VV(x) must be orthogonal to the 
generator A%x. As a consequence, the map / must be T x ^-orthogonal, due 
to the equalities 



(f(x),x)=-v{Jx,x) + V(x)\$ r =0 
(f(x), -Jx) = v X - (x, x) |§* - (Jx, VV(x)) = 0. 

Define Z„ as the subgroup of T generated by (C, C) <= x SO (2) with 
C = 2tt/ti e 50(2). Since the action of (C, C) leaves fixed the equilibrium a, 
then the isotropy group of a is the group T a x S 1 with 

Ta = <Zn- 

Thus, the orbit of a is isomorphic to the group 5*0(2). In fact, the orbit consists 
of the rotations of the equilibrium. As a consequence, the generator of the orbit 
Aia = —Ja must be in the kernel of D 2 f(a). 

3 The Liapunov-Schmidt reduction 

In order to apply the orthogonal degree of [8] , one needs to make a reduction of 
the bifurcation map to some finite space. 
The bifurcation map / has Fourier series 

iez 

where xi and gi are the Fourier modes of x and VV(x). Since —ilv(ij) is 
invertible for all big Zz/s, then one may solve xi for |2| > p and v bounded from 
below, whenever x(t) belongs to a bounded set VL in the space H 1 , that is also 
bounded uniformly in M 2 . 

In this way, the bifurcation operator / has the same zeros as the bifurcation 
function 

f(x u x 2 (xi,v),v) = ^2(~lviJxi+gi)e M , 
\l\<p 

and the linearization of the bifurcation function at some equilibrium a is 

f'(a)x 1 = {-IviJ + D 2 V{a))xie llt . 

\i\<p 
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Here X\ corresponds to the 2p+l first Fourier modes and x%(x\,v) is the result of 
applying the global implicit function theorem for functions in fi and v bounded 
from below. 

The linearization of the bifurcation function is determined by blocks M(lv) 
for I G {0, where M(y) is the matrix 

M{v) = ~vij + D 2 V{a). 

These blocks M{lv) represent the Fourier modes of the linearized equation at 
the equilibrium. 

4 Irreducible representations 

In order to apply the orthogonal degree, one needs to find the irreducible rep- 
resentation subspaces for the action of T a = Z n . 

For k G {1, n}, we define the isomorphisms Tk : C 2 — > Wk as 

T k {w) = (n- 1 / 2 e «w+J)C u , )>> . >n -i/2 e n«w+J)c u; ) with 

W k = {(e^ kI+J Kw, e n ^ m+J ^w) : w G C 2 }. 

In the paper [3], we have proved that the subspaces Wk are irreducible 
representations of the group Z„. Also, we showed that the action of (£, £) G Z„ 
on the space VFfc is given by 

p(C,0 = e ikC . 

Since the subspaces Wk are orthogonal, then the linear map 

n 
3 = 1 

is orthogonal, where w = (wi, ...,w n ). 

Since the map P rearranges the coordinates of the irreducible representa- 
tions, one has, from Schur's lemma, that 

p- l D 2 V{a)P = diag(B u B n ), 

where B>k are matrices which satisfy D 2 V(a)Tk(w) — Tk(Bkw). In the paper [3], 
we have found the blocks Bk'- they satisfy = Bk and we had the following 

result: 

Define afe and "fk as 

a/- = 4 cos £ sin 2 kQ/2 and 7^ = 2sinfc£sin£. 

Then, the blocks Bk are 

B fc = -a k I + lk {iJ) + 2[i 2 ti(\fi\ 2 )diag{l, 0). 
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For the linearization of the equation one has that 



P 1 M{v)P — diag(mx(y), m n (v)). 



Thus, we find the matrices m^iy) in terms of the blocks Bk as 



m k (v) 



v{iJ) + B k for k e {1, n}. 



The action of (C,C<?) £ t x S 1 on W fe is p{C,C,<p) = e ik ^e a ^. Therefore, 
the isotropy group of the space W k is 



5 Bifurcation theorem 

The orthogonal degree is defined for orthogonal maps that are non-zero on the 
boundary of some open bounded invariant set. The degree is made of integers, 
one for each orbit type, and it has all the properties of the usual Brouwer degree. 
Hence, if one of the integers is non-zero, then the map has a zero corresponding 
to the orbit type of that integer. In addition, the degree is invariant under 
orthogonal deformations that are non-zero on the boundary. The degree has 
other properties such as sum, products and suspensions. For instance, the 
degree of two pieces of the set is the sum of the degrees. The interested reader 
may consult [8], p] and [13] for more details on equivariant degree and degree 
for gradient maps. 

Now, if one has an isolated orbit, then its linearization at one point of the 
orbit xq has a block diagonal structure, due to Schur's lemma, where the isotropy 
subgroup of xq acts as Z n or as S 1 . Therefore, the orthogonal index of the orbit 
is given by the signs of the determinants of the sub-matrices where the action 
is as Z„, for n = 1 and n = 2, and the Morse indices of the submatrices 
where the action is as S . In particular, for problems with a parameter, if 
the orthogonal index changes at some value of the parameter, one will have 
bifurcation of solutions with the corresponding orbit type. Here, the parameter 
is the frequency v. 

The fixed point subspace of the isotropy group T s x S 1 corresponds to the 
block m„(0) = B n . Since the generator of the kernel is A\d = T„(— n 1 /^), 
then 62 must be in the kernel of m n (0). 

Following [5], one defines a to be the sign of m„(0) in the orthogonal subspace 
to C2- Since B n = p?h! (/j, 2 )diag(l,0), then 



We have proved, in [3], that m/ c (0) = B k is invertible except for a point \x k for 
k = 1, solution of fj, 2 h'(fi 2 ) — 5k, with 



Zn(*) = <(C,C,-*0>- 



a = sgn(ei B n e\ 



) = sgn(h'(n 2 )). 



5k = (4 - 1 t)/{2a k ). 
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Definition 4 Following '81, we define 

Vk{vo) = a{n k (v - p)) - n k (v + p)}, 

where n k (v) is the Morse index of m k {v). 

This number corresponds to the jump of the orthogonal index at isq. Then, 
from the results of [5], we can state the following theorem for p ^ pi, p n -\. 

Theorem 5 If rj k {v k ) is different from zero, then the relative equilibrium has a 
global bifurcation of periodic solutions from 2ir/v k with isotropy group Z„(fc). 

The solutions with isotropy group Z„(fc) must satisfy the symmetries 

u j+1 (t) = e^ Ul (t+jkO. 

The norms of these oscillators are related by the formula rj + i(t) = r%(t + 
jkQ), in particular, for k — n, all oscillate in an identical fashion. If k divides n, 
then one will have k equal traveling waves, each one formed by n/k oscillators. 

By global bifurcation, we mean that the branch goes to infinity in norm or 
period (for these possibilities we say that the bifurcation is non admissible) or, 
if not, then the sum of the above jumps, over all the bifurcation points, is zero. 

A complete description of these solutions may be found in the paper [5] . 

6 Spectral analysis 

The m k {y) have real eigenvalues, so the matrices m n - k {v) and m k (—v) have 
the same spectrum due to the equality m n ^ k (v) = fhf.(— v)- As a consequence, 
the Morse numbers satisfy 

n n -k(v) = n k {-u). 

For n — 4, one has that ctk — and the determinant of mk(v) is — (7fe — ^) 2 , 
and there is no jump in the Morse number. Note also that, for k = n, one has 
c-n = In — and the determinant is — v 2 , that is to„(^) is invertible for v ^ 
and there is no bifurcation of periodic solutions with that symmetry. 

Since we are going to treat the case n = 3 separately, we shall suppose for 
now that n > 5. For these values of n one has that a k > and 5k < ctk/ '2. 

Proposition 6 Define v± as 

is±=j k ± yjctk (a k - 2p 2 h'{p 2 )). 

For n > 5, the matrix m k (y) changes Morse index only for v±, whenever 
p 2 h'(p 2 ) < a k /2. In this case 

f)k{v±) = ±0-; 

with a = sgn(h'(p 2 )). Furthermore, the values v± are positive only in the 
following cases: 
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(a) The value v + is positive for k e {1, — 1} and /j, 2 h'([i 2 ) < 8 k . 

(b) The values i> + and V- are positive for k € {1, .., [n/2]} and 5 k < ^ 2 h'(p?) < 

afc/2. 

Proof. Since the determinant of m k is 

d k (v) = detrrifc = -2ct k \i 2 ti (\i 2 ) + ct\ - (j k - v) 2 , 

then d k (v) is zero only for v±, if fi 2 h'(fi 2 ) < a k /2. Furthermore, the trace of 
mkiy) is 

T fe = 2fi 2 h'(\^\ 2 ) - 2a k < -a k < 0. 

Since d k {v) is positive in z^ + ) and negative in the complement, then n k (v) = 
2 in and rife(z/) = 1 in the complement. Thus, r]k{^-) = f(l — 2) and 

%(!/+)= (7(2 - 1). 

Since d k {v) is a polynomial of degree two, then is positive and V- is 
negative, if d k (0) is positive. We conclude the result for (a), that is d k (0) is 
positive for /i 2 ft/(/i 2 ) < 5 k . Furthermore, d k (0) is negative for S k < /i 2 /i'(/i 2 ), 
thus, v + and V- have the same sign. Since v± — j k ± -y/*, then t/± have the sign 
of 7^. We conclude the result for (b), that is "f k is positive for k £ {1, .., [n/2]} 
and that 7„_fe = —7*5. ■ 

Thus, for the lattice with a general potential for n > 5 we have: 

Theorem 7 For eac/i fc e {1, ...,n — 1} such that /j, 2 h'(fj, 2 ) < 5 k , the relative 
equilibrium has a global bifurcation of periodic solutions starting from the period 
2tt/v + with symmetries Z„(fc) . Furthermore, this bifurcation is non-admissible 
or goes to another equilibrium. 

For each k G {1, [n/2]} such that S k < /i 2 h'(fi 2 ) < a k /2, the equilibrium 
has two global bifurcations of periodic solutions starting from the periods 2ir/v+ 
and 2Ttjv- with symmetries Z„(fc). 

Remark 8 We have proven that the determinant of the block m k {v) is zero 
at two values for k e {l,...,n — 1} if fi 2 h'(/i 2 ) < a k /2. Thus, the block m k 
corresponds to stable solutions of the linear equation for /x 2 ft/(/i 2 ) < a k /2. On 
the other hand, detm n (v) has a double zero v = 0, but the block m n gives stable 
solutions due to the spatial symmetries of the problem. Using the fact that a k 
are increasing for k e {1, [n/2]}, we conclude that the solution a is linearly 
stable for 

H 2 ti{n 2 ) < ax/2. 

The only remaining case is n = 3. The proofs arc similar to the previous 
case, by taking the reverse inequalities. 

Proposition 9 For n = 3, the determinant, detm k , is zero only at v± for 
a k /2 < fj, 2 h'(fi 2 ). In this case one has that rj k (v±) = =pa with a = sgn(h'(n 2 )). 
Furthermore, the values v± are positive only in the following cases: 

(a) // v + is positive for k G {1, 2} and < ^i 2 h'(^, 2 ). 
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(b) If v + and v_ are positive for ct\/2 < ii 2 h'(p?) < 0. 

Proof. Since a k < 0, then detm k is zero at v± if a k /2 < /j, 2 h'(ii 2 ). In this case 
the trace of m k is 

Tfc(A) = 2 / u 2 /i'(| / u| 2 ) - 2a fc > -a fc > 0. 

Hence, n k {y) = at (zv_, and ni(z^) = 1 on the complement. Thus, rik{v-) = 
(7(1-0) and %(!/+) =<r(0-l). 

Furthermore, ^+ is positive when dfe(0) is positive. Since ct\ = —71, then 
5i = 62 = 0. We conclude the result for (a) from the fact that (4(0) is positive 
for = S k < fi 2 h'(fi 2 ) . Also, (4(0) is negative for fi 2 h'(^i 2 ) < S k — 0, then v + 
and V- have the same sign. Since v± = j k ± then i/± have the sign of j k - 
We get the result for (b), from 71 = 2 sin 2 ( > and that 72 = —71. ■ 

Hence, for the lattice and n = 3, the relative equilibrium has, for each 
k = 1,2, if < fj, h'(fjr), a global bifurcation of periodic solutions starting 
from the period 2k/v + with symmetries Z 3 (fc). Furthermore, this bifurcation is 
inadmissible or goes to another equilibrium. 

For n = 3 the relative equilibrium, for ol\/2 < fjrh'(fi ) < 0, has global 
bifurcations of periodic solutions starting from the periods 2tt/u + and 2it/i>_ 
with symmetries Zs(l). 

For n = 3, the relative equilibrium is linearly stable for a\/2 < fjrh'(fj, ). 

6.0.1 Schrodinger potential 

For the lattice for the cubic Schrodinger potential, we have h(fj.) = yU with 

h'{p, 2 ) = 1 and a = 1. 

For n > 5, one has the following cases: 

For k € {1, 2,n— 2,n— 1} we have 8 k < 0, then condition (a) is never satisfied 
and, for [i e (0, y^ajfe/2), condition (b) holds. 

For k € {3, ...,n— 3} we have 5fc > 0, then, for fi £ (0, \f8~k), condition (a) is 
satisfied, while, for fi £ (\/3fe, ^a k /2), condition (b) holds. 

Theorem 10 Thus, for the lattice with a Schrodinger potential, for n > 6, the 
relative equilibrium has, for each k £ {3, ...,n — 3} and /1 £ (0,\/~8~k), a global 
bifurcation of periodic solutions, starting from the period 2ir/v+ with symme- 
tries Z„(fc). Furthermore, this bifurcation is inadmissible or goes to another 
equilibrium. 

For n > 5 the relative equilibrium has, for each k £ {1,2} such that [i £ 
(0, \Jo.k/2) and, for each k £ {3, [n/2]} such that fi £ (y/$k, yJa k /2), two 
global bifurcations of periodic solutions starting from the periods 2t:/v + and 
2it/v- with symmetries Z n (fc). 

For n > 5, the relative equilibrium is linearly stable for the amplitudes 

fl£ (0, 
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For n = 3 and k G {1,2}, we have that condition (a) is always satisfied. For 
n = 3, the relative equilibrium has, for each k G {1,2} such that /i G (0, oo), 
a global bifurcation of periodic solutions starting from the period 2tt/i> + with 
symmetries ^s(k). Furthermore, for n — 3, this relative equilibrium is linearly 
stable for the amplitudes fi G (0,oo). 

6.0.2 Saturable potential 

For the lattice with the saturable potential, one has h{x) 
h'(fi 2 ) = -(1 + /I 2 )" 2 and a = -1. 

For n > 5 we have the following cases: 

For k G {2, ...,n — 2} one has 5k > 0, then, for each fi G (0, oo) one gets 
/j 2 h'(^i 2 ) < 5k, that is condition (a). 

For n G {16,17,...}, let G (0,1) and fi + G (l,oo) be the solutions of 
/z 2 ft/(/i 2 ) = Si. Then, for \i G (0, /x_) U (/z+,oo), one has 5i < fx 2 h'(p 2 ), that 
is condition (b) and, for G (/i_,/i + ) one has /j, 2 h'(^i 2 ) < 5\, that is condition 
(a). Furthermore, for n = {5, 15} one has the same result as before provided 
one defines and to be zero. 

Theorem 11 Thus, for the lattice with the saturable potential, for n > 5, the 
relative equilibrium has, for each fc € {2, n — 2} and /i£ (0, oo) and, for each 
k € {1, ra— 1} anrf ^ G /i+), a global bifurcation of periodic solutions starting 
from the period 2-k/v + with symmetries Z„(fc). T/iis bifurcation is inadmissible 
or goes to another equilibrium. 

Moreover, the relative equilibrium has, for each fj, G (0,/i_) U (/i+,oo) 7 £wo 
g/o&aZ bifurcations of periodic solutions starting from the periods 2t:/v + and 
2-k/v- with symmetries Z„(l). 

For n > 5, the relative equilibrium is linearly stable for all amplitudes fi G 
(0,oo). 

For n — 3, we have that ai/2 = -3/4, hence, ai/2 < -1/4 < p?h\p?) for 
all G (0, oo). Thus, condition (b) is always satisfied. For n = 3, t/ie relative 
equilibrium has, for each [i G (0, oo), global bifurcations of periodic solutions 
starting from the periods 2tt/v + and 2-k/v- with symmetries Z 3 (l). Further- 
more, for n = 3, the relative equilibrium is linearly stable for the amplitudes 
fj, G (0, oo). 
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